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Abstract 

The chemical equilibration of the parton distribution in collisions of two 
heavy nuclei at the CERN Large Hadron Collider is investigated. Initial 
conditions are obtained from a self-screened parton cascade calculation. The 
onset of transverse expansion of the system is found to impede the chemical 
equilibration. The system initially approaches chemical equilibrium, but then 
is driven away from it, when the transverse velocity becomes large. 
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It is expected that very highly energetic heavy ion collisions lead to the formation of a 
new phase of strongly interacting matter, called the quark-gluon plasma (QGP). At collider 
energies it seems reasonable to visualize the two nuclei as two clouds of valence and sea 
partons which pass through each other and interact JTJ. The multiple parton collisions are 
thought to produce a dense plasma of quarks and gluons. After its formation this plasma 
will expand, cool, and become more dilute. If quantum chromodynamics admits a first- 
order deconfinement or chiral symmetry phase transition, it is likely that the system will 
pass through a mixed phase of quarks, gluons, and hadrons, before the hadrons lose thermal 
contact and stream freely towards the detectors. 

Several questions about the structure of the matter formed in these nuclear collisions 
arise. Does the initial partonic system attain kinetic equilibrium? Probably yes, as the 
initial parton density is large, forcing the partons to suffer many collisions in a very short 
time || . Does it attain chemical equilibration? This will depend on the time available [0-[5| 
to the partonic system before it converts into a mixed phase. The time available for this is 
perhaps too short (3-5 fm/c) at the energies (-^/s < 100 GeV/nucleon) to be reached at the 
Relativistic Heavy Ion Collider (RHIC). At the energies (y/s < 3 TeV/nucleon) that will be 
achieved at the CERN Large Hadron Collider (LHC) this time could be large (more than 
10 fm/c). If one considered only a longitudinal expansion of the system, the QGP formed 
at LHC energies would approach chemical equilibrium very closely, due to the higher initial 
temperature predicted to be attained there. However, the life-time of the plasma would also 
be large enough to allow a rarefaction wave from the surface of the plasma to approach the 
center. We shall see that this could have interesting consequences for the evolution of the 
plasma. 

It has recently been shown || that color screening provides a mechanism for the elimina- 
tion of infrared divergences in the partonic cascades following the interaction of two heavy 
nuclei. Early, hard scatterings produce a medium which screens the longer ranged color 
fields associated with softer interactions. When two heavy nuclei collide at high enough 
energy, the screening occurs on a length scale where perturbative QCD still applies. This 
approach yields predictions for the initial conditions of the forming QGP without the need 
for cut-off parameters. The resulting values to be adopted in the present study are listed 
in Table I. They show that the QGP initially is extremely hot, with temperatures around 1 
GeV, but not chemically equilibrated. 

Once the kinetic equilibrium has been attained we assume that the system can be de- 
scribed by the equation for conservation of energy-momentum of an ideal fluid: 

dfjT^ = , = (e + P)u»u v + Pg^ (1) 

where e is the energy density and P is the pressure measured in the frame comoving with the 
fluid. We use natural units (h — c — 1). The four- velocity vector for the fluid satisfies the 
constraint u 2 = —1. For a partially equilibrated plasma of massless particles the equation 
of state can be written as 

e = 3P=[a 2 \ g + b 2 (\ q + \ q )}T* (2) 

where a 2 = 87r 2 /15, b 2 = 77r 2 iVj/40, Nf « 2.5 is the number of dynamical quark flavors, 
and Afc is the fugacity for the parton species k. We take the speed of sound c s = l/y/3 as 
implied by @ and solve the hydrodynamic equations ([!]) for a system undergoing a boost 
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invariant longitudinal (z-axis) and cylindrically symmetric transverse expansion (see Ref. 
for details) with initial conditions obtained earlier. It is sufficient to solve the problem 
for z = 0, because of the assumption of boost invariance. Confining ourselves || to the 
dominant reactions gg «-> ggg, and gg <-> qq for the equilibration of the partons, we write 
the master equations governing the parton densities as || 

dpingvP) = n g (R 2 ^ 3 ~ - {ngRg-^ q - n q R q ^ g ) 

dfj,(n q u^) = d^riqU^) = rigRg^g - n q R q ^ g , (3) 

in an obvious notation. These equations lead to partial differential equations for the fugac- 
ities: 

Y g dAg + ^d r \ 9 + ^( 7 r 3 ) + ^( 7 t 3 ) 

+ "/d r v r + 7 + 

= R 3 (l - X g ) - 2R 2 U-^j , 



(4) 




where R 2 and R3 are related to the rates appearing in (|3]) and are given by, 

R 2 w 0.24iV / a^A g Tln(1.65/Q; s A g ) 

R 3 = 1.2a 2 s T(2X 9 - X 2 ) 1/2 , (5) 



v r is the transverse velocity and 7 = — vf. These will reduce to (34, 35) of Biro et 

al. 0, if there is no transverse expansion of the system. The hydrodynamic equation ([I]) is 
numerically solved to get e(r, t) and v r (r,t) which are then used to solve the equations for 
the fugacities 0. We have verified that our results near r = closely follow the results for 
a purely longitudinal expansion, till the time when the rarefaction wave has not reached the 
center. We have solved the hydrodynamic equations by assuming that the initial transverse 
velocities are zero. 

Recall that for a boost invariant longitudinal expansion ([I]) provides that er 4 / 3 is a 
constant |§. Thus in order to bring out the consequences of transverse expansion in a 
transparent manner, we plot the constant energy density contours e(r, t) = ei/N 4 ^ 3 for 
N=l,6,ll..,71 (Figure la), which would be lines parallel to the N=l case extending upto 
r = Rt and separated by AN t{ , if there were no transverse expansion of the system. In the 
present case N=71 corresponds to an energy density [ 10 1 of 1.45 GeV/fm 3 below which we 
assume the description in terms of a perturbative QCD plasma to fail. We have also shown a 
line r = Rt — c s t which indicates the radial distance of the region affected by the rarefaction 
wave at any given time t. We see that almost the entire fluid is likely to be affected by the 
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flow as the life time of the system is sufficiently large. A more rapid cooling of the system 
is also indicated by the decreasing separation of the contours at later times at, say, r = 0. 
This will affect the evolution of the chemical equilibration in several competing ways. 

Thus in Figs, lb and lc we show the fugacities for the gluons and the quarks, respectively, 
along some of the above constant energy density contours. Several interesting observations 
can be made. We note that the quark fugacities lag behind the gluon fugacities at all times 
and all radial distances. This is not surprising, considering the lower starting values for \ q . 
Recalling that increasing N denotes passage of time, we see that the fugacities at moderate 
values of r first increase with passage of time. However beyond N = 31, (t ~ 8 fm/c, near 
r = 0) when the transverse expansion of the system is large, it starts decreasing again, except 
at r = where it continues to increase, as all the radial derivatives in (f|) are negligible there. 

This can be understood by noting that the evolution of the partonic density which is 
governed by d^nu 11 ) consists of two terms: u^d^n which is the rate of change of the density 
in the comoving frame, and nd^u^ which is the rate of change of the density due to the 
expansion of the fluid element in the comoving frame ||. Once the transverse expansion of 
the fluid starts developing, the second term and also the radial derivatives in (^) grow very 
rapidly and drive the system away from chemical equilibrium. 

The fugacities decrease monotonically with increase in r which is not surprising because 
the transverse flow effects grow with increasing r and also the time available to the system 
in the QGP phase shrinks at larger r. This can have very interesting consequences, the 
most intriguing of which is the possibility of the formation of a hot and rare partonic system 
near the surface as compared to a cooler but denser partonic system at the center. Thus, 
e.g., we estimate that the temperature at the end of our QGP phase is about 200 MeV at 
r=6 fm, but only about 160 MeV at r=0 fm. This variation in the temperature gets more 
severe at RHIC energies [IT|, as the life-time of the QGP phase is much smaller there, thus 
terminating the journey of the partonic system far away from chemical equilibrium. One 
may create arbitrary competitions between the life-time and the transverse flow effects by 
considering systems of varying dimensions and initial conditions. 

It is not clear whether such a system will go through a mixed phase, say by nucleation of 
a hadronic droplet |12[] , at all radial distances even if QCD admits a first order quark-hadron 
phase transition. Lowering the energy density at which we terminate our calculations will 
only sharpen this difference as the transverse flow will further enhance the gradients. 

The development of the transverse flow and its consequences on the chemical equilibration 
are sensitively dependent on the evolution of the transverse velocity of the fluid v r (r,t). In 
Fig. 2 we have shown the variation of the transverse velocity along some of the constant 
energy density contours shown in Fig. la. We see that the transverse velocity develops 
rapidly and spreads to smaller radii very quickly, as the speed of sound which drives the 
pressure is assumed to be rather large. It is conceivable that the speed of sound is reduced 



when the energy density gets smaller JT3[. The life-time of the system would, however, 
be larger then, and thus the system remains prone to the consequences of the transverse 
expansion, as seen earlier. 

If we now assume that the matter goes to a mixed phase (at small radii, at least), the 
transverse velocity attained at the end of the QGP phase will stay constant for a compar- 
atively long time, and it may be possible to devise means to actually determine it from 
the measured particle spectra. Since D-mesons and even B-mesons are expected to be co- 
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piously produced at LHC energies, the decay leptons from these heavy mesons could serve 
as indicators of the collective transverse flow, because the leptons will inherit a substantial 
fraction of the significant kinetic energy carried by their parents. Due to their large mass 
and relatively small cross sections, the heavy mesons cannot be accelerated to acquire such 
high velocities in the hadronic phase. 



The transverse flow will also affect the spectra of photons and dileptons mM. We find 



II | that the yield of single photons with large transverse momenta increases by a factor of 



almost 5 due to the transverse expansion, as compared to the situation when there is only 
a longitudinal expansion. Similar trends are also seen in the transverse mass distribution 
of dileptons. We have seen that gluon fugacities are larger than the quark fugacities. This 
provides that the annihilation processes (oc X q Xq), give only a small contribution as compared 
to the Compton processes (oc X g X q ) to the production of single photons, and can be neglected. 
Dileptons on the other hand are produced from quark- ant iquark annihilation. A comparison 
of these two measurements could be a valuable source of information on the quark and gluon 



fugacities of the system [15 



An experimental measurement of the strong transverse flow predicted here would provide 
two important pieces of evidence about the properties of the dense matter produced in 
ultrarelativistic heavy ion reactions: that the matter is formed at very high initial pressure, 
and that this pressure is maintained for a substantial period of time. If the "rarified" 
partonic matter hadronizes before it attains chemical equilibrium, as is predicted at RHIC 
energies, it is possible that a mixed phase or a interacting hadronic phase is not created. The 
emitted hadrons would then reflect the thermal characteristics of the quark phase from which 
they emerge. We are currently investigating whether another initial parton distribution 
would lead to a significantly different transverse flow pattern that could be experimentally 
discriminated from the one predicted here. 
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TABLES 



n = 0.25 fm/c 


RHIC 


LHC 


€t (GeV/fm 3 ) 


61.4 


425 


T % (GeV) 


0.668 


1.02 


x» 
A 9 


0.34 


0.43 


yii) 


0.064 


0.082 



TABLE I. Initial conditions for the hydrodynamical expansion phase at RHIC and LHC. The 
initial time is taken as Tj = 0.25 fm/c. 
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FIG. 1. (a) Constant energy density contours for a transversely expanding quark-gluon plasma 
likely to be created in Au+Au collisions at LHC. For a boost invariant longitudinal expansion these 
contours would be equidistant (At = 5t{ at r =0) lines parallel to N = 1. The dashed lines gives 
the distance from the axis where the rarefaction wave has arrived at time t, and beyond which the 
fluid is strongly affected by the flow, (b) Gluon fugacity along some of the constant energy density 
contours, (c) Quark fugacities along some of the constant energy density contours. 
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FIG. 2. Transverse velocities along some of the constant energy density contours. 
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